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ANALYSIS  OF  A  FLUID-LOADED  THICK  PLATE 


INTRODUCTION 

The  physics  of  a  thick  plate  with  fluid  loading  on  both  sides  provides  the  theoretical  basis  for 
insertion  loss  and  echo  reduction  tests,  both  of  which  are  typically  used  to  determine  how 
efficiently  a  material  transmits  or  reflects  energy.  Such  testing  is  conducted  by  insonifying  a 
submerged,  slab-shaped  sample  and  then  measuring  the  transmitted  and  reflected  sound  pressure. 
Based  on  the  sound  pressure  level  of  the  incident  field,  the  insertion  loss  and  echo  reduction 
quantities  can  be  calculated.  When  these  experiments  are  performed  in  a  small  tank,  the  wall 
motion  of  the  sample  is  sometimes  measured,  with  the  fluid  pressure  then  calculated  based  on 
this  displacement.  The  above  measurements  and  corresponding  theory  are  important  to  the  U.S. 
Navy  because  they  help  designers  develop  the  most  effective  acoustic  sonar  windows,  sonar  and 
ship  baffles,  and  anechoic  marine  coatings. 

Central  to  the  study  and  understanding  of  insertion  loss  and  echo  reduction  is  plate  theory, 
which  has  been  researched  extensively  for  many  years.  Thin  plate  theory,1  studied  in  1973,  is  a 
simplified  version  that  fails  to  accurately  incorporate  dynamic  response  when  the  sample  is  thick 
compared  to  a  wavelength.  In  contrast,  thick  plate  theory2  usually  incorporates  all  the  dynamics 
of  the  plate  and  is  normally  used  when  the  sample  is  on  the  order  of  a  wavelength  of  energy  in 
the  structure.  More  complex  investigations  have  analyzed  the  dispersion  curve  for  the  plate 
without  fluid  loading3'5  or  for  the  plate  in  contact  with  a  continuous  fluid  on  one  or  both  sides.6'11 
To  a  lesser  extent,  papers  have  been  published  that  examine  plate  response  to  various  other 
loading  configurations.  For  example,  studies  have  explored  the  radiation  efficiency  of  infinite 
fluid-loaded  plates  subjected  to  point  loads,12  calculated  the  corresponding  transfer  functions  for 
thin  plate  models  coupled  to  fluid  loading,13  and  determined  mode  shapes  for  a  thick  plate  with 
finite  depth  that  is  loaded  by  fluid  on  both  sides.14  During  these  investigations,  the  response  of 
thick-walled  plates  has  been  typically  left  as  an  open-form  solution  that  involves  a  matrix  inverse 
at  a  specific  wavenumber  and  frequency.  Although  closed-form  solutions  have  been  previously 
derived,  such  analyses  have  been  performed  only  at  zero  wavenumbers  and  have  not  been 
extended  to  nonzero  wavenumbers. 


1 


This  report  derives  the  equations  of  motion  of  an  infinite  thick  plate  coupled  on  one  or  both 
sides  with  fluid  loading  as  it  is  excited  with  a  continuous  forcing  function.  The  equations  of 
motion  are  formulated  into  a  four-by-four  system  of  linear  equations  with  wave  propagation 
coefficients  as  the  unknown  terms.  Once  the  system  matrix  is  known,  the  dispersion  equations 
—  derived  in  closed-form  expressions  from  the  determinant  of  the  matrix  —  explicitly  show  the 
effects  of  the  fluid  loading.  Calculated  next  are  the  closed-form  transfer  functions  of  plate 
motion  divided  by  source  excitation,  which  contain  the  plate  and  the  fluid  terms  separately. 
Based  on  these  transfer  functions,  the  displacement  shape  of  the  plate  modes  is  studied  with 
respect  to  fluid  loading  on  one  or  both  sides  of  the  plate.  Insertion  loss  and  echo  reduction  are 
then  calculated  for  the  system  at  nonzero  wavenumbers  using  closed-form  solutions  of  the 
pressure  fields  that  are  determined  based  on  plate  wall  motion,  with  these  results  compared  to  the 
previously  available  values  at  zero  wavenumber. 


SYSTEM  MODEL 

As  shown  in  figures  1  and  2,  the  system  model  is  a  thick  plate  in  contact  on  one  or  both  sides 
with  a  fluid  that  exerts  a  continuous  excitation  pressure  on  the  plate.  The  model  configurations, 
referred  to  as  single  and  double  fluid-loaded  plates,  are  based  on  the  following  assumptions: 

(1)  the  forcing  function  acting  on  the  plate  is  a  plane  wave  with  definite  wavenumber  and 
frequency  content,  (2)  the  corresponding  response  of  the  plate  is  at  a  definite  wavenumber  and 
frequency,  (3)  motion  is  normal  and  tangential  to  the  plate  in  one  direction  (two-dimensional 
system),  (4)  the  plate  has  an  infinite  spatial  extent,  (5)  the  particle  motion  and  pressure  response 
is  linear,  and  (6)  the  fluid  medium  has  no  loss.  For  the  case  where  the  fluid  is  on  both  sides  of 
the  plate,  each  fluid  has  the  same  acoustic  properties. 

For  the  single  and  double  fluid-loaded  plates,  the  acoustic  pressure  in  the  fluid  on  the 
excitation  side  of  the  plate  is  governed  by  the  wave  equation  and  is  written  in  Cartesian 
coordinates  as15 


#2P\(x,z,t) 

dz2 


d2P\(x>z,t)  1  ^2p,(x,z,0 

dx2  c)  dt2 


=  0  , 


0) 


2 


Figure  1.  Single  Fluid-Loaded  Plate 


Fluid  1 


Fluid  2 


Figure  2.  Double  Fluid-Loaded  Plate 

where  px{x,z,t)  is  the  pressure  (N/m2),  with  the  subscript  1  denoting  the  fluid  on  the  excitation 
side  of  the  plate;  z  is  the  spatial  location  (m)  normal  to  the  plate;  jc  is  the  spatial  location  (m) 
tangential  to  the  plate;  cf  is  the  compressional  wave  speed  of  the  fluid  (m/s);  and  t  is  time  (s). 
The  coordinate  system  of  this  configuration  is  shown  in  figure  3.  Note  that  the  use  of  this 
orientation  results  in  b  =  0  and  a  having  a  value  less  than  zero.  Furthermore,  the  thickness  of  the 
plate,  h,  is  a  positive  value,  and  the  motion  of  the  plate  is  governed  by  the  equation16 

//V2«  +  (/l  +  ^)W.u  =/>0,  (2) 

where  p  is  the  density  (kg/m3),  X  and  p  are  the  Lame  constants  (N/m2),  •  denotes  a  vector  dot 
product,  and  u  is  the  Cartesian  coordinate  displacement  vector  of  the  plate.  For  the  double 


Figure  3.  Coordinate  System  of  Model 


fluid-loaded  plate,  the  acoustic  pressure  in  the  fluid  opposite  the  excitation  side  of  the  plate  is 
governed  by  the  wave  equation  and  is  written  in  Cartesian  coordinates  as 

d\p2(x,z,t)  02p2(x,z,t)  1  d2  p2(x,  z,  t)  n 

+— ^  — =  0’  <3> 

where  p2(x,z,t )  is  the  pressure  (N/m2)  and  the  subscript  2  denotes  the  fluid  opposite  the 
excitation  side  of  the  plate.  For  the  single  and  double  fluid-loaded  plates,  the  interface  between 
the  first  fluid  and  the  surface  of  the  plate  at  z  =  b  satisfies  the  linear  momentum  equation,  which 


„  d2uz(x,b,t)  dpfx,b,t) 

Pf  — W2  = - (4) 

where  pf  is  the  density  of  the  fluid  (kg/m3).  For  the  double  fluid-loaded  plate,  the  interface 
between  the  second  fluid  and  the  surface  of  the  plate  at  z  =  a  also  satisfies  the  linear  momentum 
equation  and  is  written  as 

„  dp2{x,a,t) 

p‘  —J?  = — <5> 

B)uations  ( I )  through  (5)  are  the  governing  partial  differential  equations  of  the  single  and  double 
fluid-loaded  plate  systems. 
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Equations  (1)  through  (3)  are  now  transformed  from  partial  differential  equations  into 
ordinary  differential  equations  and  then  into  algebraic  expressions.  The  acoustic  pressure  in 
equation  (1)  is  modeled  as  a  function  at  definite  wavenumber  and  frequency  by 

(x,  z,t)  =  P]  (z,  kx ,  co)  exp(i^x)  exp(i&>?)  ,  (6) 

where  co  is  frequency  (rad/s),  kx  is  the  spatial  wavenumber  in  the  x-direction  (rad/m),  and  i  is  the 
square  root  of-1.  If  the  pressure  in  the  fluid  is  generated  by  an  acoustic  plane  wave,  the  spatial 
wavenumber  is  given  by 

K=—  sin(0),  (7) 

cf 

where  6  is  the  angle  of  incidence  (rad)  of  the  incoming  acoustic  wave,  with  6=  0  corresponding 
to  excitation  normal  to  the  plate  (or  broadside  excitation).  Wavenumbers  larger  than  co/cf  are 
possible  and  are  typically  generated  from  turbulent  fluid  loading  or  structural  wave  loading. 
Inserting  equation  (6)  into  equation  (1)  and  solving  the  resulting  ordinary  differential  equation 
yields 

Pl ( z , kx,0))  =  H(kx , a) exp(iyz)  +  Pe (co) exp(-i^z)  .  (8) 

In  equation  (8),  the  first  term  on  the  right-hand  side  represents  the  reflected  (or  reradiated) 
pressure  field  and  the  second  term  represents  the  applied  incident  pressure  field  (the  excitation  or 
forcing  function)  acting  on  the  plate.  To  be  exact,  the  term  H(kx,co)  is  the  wave  propagation 
coefficient  of  the  reflected  pressure  field,  and  the  term  Pe(co)  is  the  excitation  (or  source)  level. 
Furthermore, 


r  = 


UY 


\cf  J 


where  y  is  the  wavenumber  of  the  acoustic  pressure  in  the  fluid  (rad/m). 


(9) 
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Equation  (2)  is  manipulated  by  writing  the  Cartesian  coordinate  displacement  vector  u  as 
ux(x,y,z,t) 

u  =  <uy(x,y,z,t)  >  ,  (io) 

u:(x,y,z,t) 

withy  denoting  the  direction  into  the  plate.  The  symbol  V,  from  equation  (2),  is  the  gradient 
vector  differential  operator  written  in  three-dimensional  Cartesian  coordinates  as18 


„  d  .  d  .  d  . 

V  =—ix  +—  /„  +  —  l.  , 

dx  dy  }  dz 


(11) 


with  ix  denoting  the  unit  vector  in  the  x-direction,  iy  denoting  the  unit  vector  in  the  y-direction, 
and  4  denoting  the  unit  vector  in  the  z-direction.  The  symbol  V2  is  the  three-dimensional 
Laplace  operator  operating  on  vector  u  as 


V2u  =  V2uxix  +  V2uyiy  +  V2uJz 


the  term  V  •  u  is  the  divergence,  which  is  equal  to 


dux  du 

V  •  u  =  — -  + L  + 

dx  dy 


du^ 

dz 


(14) 


The  displacement  vector  u  is  written  as 
u  =  V^+V  x  ip  , 


(15) 


where  (j>  is  a  dilatational  scalar  potential,  x  denotes  a  vector  crossproduct,  and  p  is  an 
equivoluminal  vector  potential  expressed  as 
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(16) 


¥  = 


y/x(x,y,z,t) 
y/y(x,y,z,t ) 
Yz(x,y,z,t) 


The  formulation  is  a  two-dimensional  problem;  thus,  y=  0  and  d(-)  /  dy  =  0 .  Expanding 
equation  (15)  and  breaking  the  displacement  vector  into  its  individual  nonzero  terms  yields 


dx  dz 


(17) 


and 


uz(x,z,t) 


(fyjXjZ, t)  |  ^(*,z,Q 
dz  dx 


(18) 


Equations  (17)  and  (18)  are  next  inserted  into  equation  (2),  which  results  in 


c2dV2<f>(x,z,t) 


d2<f)(x,z,t ) 
~dt2 


(19) 


and 


c2sV2¥y(x,z,t) 


d2y/y{x,z,t) 

~d? 


(20) 


where  equation  (19)  corresponds  to  the  dilatational  component  and  equation  (20)  corresponds  to 
the  shear  component  of  the  displacement  field.19  Correspondingly,  the  constants  and  cs  are 
the  complex  dilatational  and  shear  wave  speeds,  respectively,  and  are  determined  by 


I A  +  2ju 

V  P 


(21) 


(22) 
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The  relationships  of  the  Lame  constants  to  the  Young’s  (compressional)  and  shear  moduli  are 
shown  as 

.  Ev 

~  (l  +  u)(l-2w)  (23) 


M  =  G  = 


2(1  +  v)  ’ 


where  E  is  the  Young’s  modulus  (N/m2),  G  is  the  shear  modulus  (N/m2),  and  u  is  the  Poisson’ 
ratio  of  the  material  (dimensionless). 

The  conditions  of  infinite  length  and  steady-state  response  are  now  imposed,  allowing  the 
scalar  and  vector  potential  to  be  written  as 

(j){x,z,t)  =  O(z)  Qxp(ikxx)  exp(i  co  t)  (25) 


y/y(x,z,t)  =  xP(z)exp(i^x)exp(i®/) 


Inserting  equation  (25)  into  equation  (19)  yields 


d 20(z)  , 


+  a  <t>(z)  =  0 


where 


M-kl  , 


Kd  ~ 


Inserting  equation  (26)  into  equation  (20)  produces 
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+  /?2vF(z)  =  0  , 


(30) 


d2x¥(z) 

dz 2 


where 

with 


(31) 


(32) 


The  solution  to  equation  (27)  is 

Q>{z)  =  A(kx,co) exp (i a z)  +  B{kx,co) exp (-iaz)  ,  (33) 

and  the  solution  to  equation  (30)  is 

'P(z)  =  C(kx ,  co)cxp(ij3z)  +  D(kx ,  a>)exp(-i/3z)  ,  (34) 

where  A(kx ,  co) ,  B(kx ,  co) ,  C(kx ,  co) ,  and  D(kx ,  co)  are  wave  propagation  constants  of  the  plate 
and  are  determined  below.  The  displacements,  now  written  as  functions  of  the  unknown 
constants  using  the  expressions  in  equations  (17)  and  (18),  are 

ux  (x,  z,  t )  =  Ux  (kx ,  z,  tf))exp(i£xx)exp(iflrf) 

=  \A{kx ,  co)ikx  exp(i  az)  +  B(kx ,  o)ikx  exp(-iaz) 

-  C(kx ,  co) i/3 exp(i/?z)  +  D(kx ,  co)\(3 exp(-i/?z)]exp(i£xx)  exp(i<wt)  (35) 

and 


uz  (x,  z,  t)  =  Uz  (kx ,  z,  co)Qxp(\kxx)txp(icot) 

=  \A(kx ,  <y)i  a  exp(i  az)  -  B(kx ,  co)  i  a  exp(-i  az) 

+  C(kx ,  co)ikx  exp(i fiz)  +  D(kx ,  co)\kx  exp(-i/? z)]exp(iArxx) exp(iot) .  (35) 
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The  normal  stress  at  the  top  of  the  plate  (z  =  b)  is  equal  to  the  opposite  of  the  pressure  in  the 
fluid  and  is  expressed  as 


(*,  b,f)  =  (X  +  2  ft)  — ■ (*’b,,)  +  X  — = 

dz  dx 


The  tangential  stress  at  the  top  of  the  plate  is  zero  and  is  written  as 


dz  dx 


For  the  double  fluid-loaded  plate,  the  normal  stress  at  the  bottom  of  the  plate  (z  =  a)  is  equal  to 
the  opposite  of  the  pressure  in  the  fluid.  This  expression  is 

«•■■(*.«,<)  =  (X  *  2 _  (39) 

where  p2(x,a,t)  represents  the  transmitted  (or  radiated)  acoustic  pressure  in  the  fluid  field  on 
the  opposite  side  of  the  acoustic  excitation.  For  the  single  fluid-loaded  plate,  p2(x,z,t )  =  0  in 
equation  (39).  The  tangential  stress  at  the  bottom  of  the  plate  is  zero,  with  this  equation  written 


dz  dx 


For  the  double  fluid-loaded  plate,  the  acoustic  pressure  in  equation  (3)  is  modeled  as  a 
function  at  definite  wavenumber  and  frequency,  resulting  in 

=  P2{z,kx,(o)exp(\kxx)cxp(\(ot)  .  (41) 

Inserting  equation  (41)  into  equation  (3)  and  solving  the  resulting  ordinary  differential  equation 
yields 


P2 (z, kx,co)  =  K (kx , co) exp(-iyz)  , 
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which  is  the  transmitted  (or  outgoing)  acoustic  energy  in  the  second  fluid.  The  term  K(kx,  co)  is 
the  wave  propagation  coefficient  of  the  transmitted  pressure  field.  Note  that  there  is  no  incoming 
wave  energy  on  this  side  of  the  plate  and  thus  only  one  exponential  term  is  present. 

Assembling  equations  (37)  through  (40);  incorporating  equations  (4),  (5),  (8),  and  (42);  and 
letting  b  =  0  yields  the  four-by-four  system  of  linear  equations  that  model  the  system  as  follows: 

Ax  =  b  ,  (43) 

where  the  entries  of  equation  (43)  are 


AUp  =  -a2A-2a2/j-Jlk 2  , 

(44) 

_pfco2a 

^lli  _  A\\d  ~  5 

y 

(45) 

> 

ii 

*5* 

+ 

> 

(46) 

> 

to 

II 

> 

1 

> 

ST 

(47) 

AaP  =  2*Jfi  , 

(48) 

,  A  Pf<*2K 

A13s  ~  A13 d  ~  > 

r 

(49) 

A\3  ~  ~A\3 p  +  A13 s  ’ 

(50) 

A14  =  A\3 p  +  A\ 3s  ’ 

(51) 

A2l  =  -2 /Axa  , 

(52) 

A  22  =  ~A21  ’ 

(53) 

A23=pP2  , 

(54) 

A  -A 

yi24  —  si23  9 

(55) 
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Aii=(Aup-And)c\p(\aa)  , 

(56) 

^32=UiP  +  ^i  Jexp(-iaa)  , 

(57) 

^33  =  (-  AnP  ~  And)exp(ij3a)  , 

(58) 

A 34  =  (aUp  -  A13d)exp(-ij3a)  , 

(59) 

A4,  =  A2]  exp(iaa)  , 

(60) 

A42  =  -A2lexp(-iaa)  , 

(61) 

A4 3  =  A2]  exp(i fid)  , 

(62) 

A  44  =  A2J  exp(-i pa)  , 

(63) 

*n  =A(kx,a ))  , 

(64) 

x2l=B{kx,co)  , 

(65) 

*3i  =C(kx,a))  , 

(66) 

*4i  =  D(kx,co)  , 

(67) 

b\\p  =  , 

(68) 

bUs=-Pe(M)  , 

(69) 

bU  ~  bllp  +  bUs  > 

(70) 

o 

II 

(71) 

631  -  0  , 

(72) 

o 

II 

(73) 

In  equations  (44)  through  (73),  the  subscript  p  corresponds  to  terms  related  to  the  plate,  the 
subscript  s  corresponds  to  the  case  of  both  the  single  and  double  fluid-loaded  plates,  and  the 
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subscript  d  corresponds  only  to  the  case  of  the  double  fluid-loaded  plate.  To  model  the  behavior 
of  the  thick  plate  without  the  fluid  loads,  the  terms  with  the  subscript  s  and  d  are  set  equal  to 
zero.  To  model  the  behavior  of  the  fluid  plate  with  a  single  fluid  load,  the  terms  with  the 
subscript  d  are  set  equal  to  zero. 

For  the  single  and  double  fluid-loaded  plate,  the  reflected  acoustic  field  on  the  excitation  side 
of  the  plate  is 


PR(kx,zb,co) 


Uz(kx,b,G>)  + 1 


exp(i yzb)  , 


(74) 


with  zb  being  the  position  where  the  field  is  evaluated  (m).  The  total  pressure  field  on  the 
excitation  side  is  a  sum  of  the  reflected  field  and  the  spatially  phase-shifted  excitation  level.  It  is 
written  as 


PTo,ai (K ,zb,<»)  =  PR (kx ,zb,(d)  +  Pe O) exp(-i yzb )  .  (75) 

For  the  double  fluid-loaded  plate,  the  transmitted  pressure  field  on  the  opposite  side  of  the 
excitation  is 


PT(kx,za,co) 


U2(kx,a,0) ) 


exp(-i/zj  , 


(76) 


with  za  being  the  position  where  the  field  is  evaluated  (m).  The  insertion  loss  (IL)  of  the  double 
fluid-loaded  plate  is  then  calculated  using 


IL(kx,a>)  =  201og10 


PM) 

pAK,co)  ’ 


(77) 


where  IL(kx,co )  is  in  units  of  decibels.  The  echo  reduction  (ER)  is  calculated  using20 
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ER{kx,co)  =  20  log, 


PM) 

PR(kxM 


where  ER(kx,o>)  is  in  units  of  decibels. 


DISPERSION  EQUATIONS 

The  dispersion  equation  is  an  equation  whose  zeros  correspond  to  single-mode  propagation 
in  the  structure.  This  function  is  proportional  to  the  determinant  of  A  in  equation  (43).  For  the 
case  of  the  single  fluid-loaded  plate,  the  equation  is  written  as 

&s(k. x’*0)  -  P\ (kx > cos(a h) cos(ph)  +  fl(kx,eo)cos(ah)sin(j3h) 

+  f2{kx,co)sm{ah)cos{ph)  +  p2(kx,co)s\n{ah)sm{ph)- pl(kx,of)  ,  (79) 


where 


Pi (kx , at)  =  -8 apk\ (, p 2  -  k2x  )2  , 


Mkz,co)  =  ipf(rpy  a(P2  ~ kx)2 (p2  +  k2x)2  , 


f2(kx,<0)  =  4ipf(rpy'a2pk2(p2  +  k2)2  , 


Pi(kx,co)  =  (P2  -k2)4  +\6a2p2k4x  . 


For  the  case  of  the  double  fluid-loaded  plate,  the  equation  is  written  as 


Ad(kx,co)  =  Pl  (kx ,  a)  cos(ah)  cos  (ph)  +  2  fx  (kx ,  a)  cos  (ah)  sin  {ph) 
+  2/2  (kx ,  a)  sin(a  h)  cos  (ph) 

+  Ip  2  (k* > +  A  ( kx » ®)]sin(«/2)  sin  {ph)  -  px  (kx ,  co)  , 
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where 


MKM  =  p2f(YPY2a\P2  +  klY  •  (85) 

In  equations  (79)  and  (84),  the  p  constants  correspond  to  the  plate  and  the / constants  correspond 
to  the  fluid  loads.  In  the  absence  of  fluid  loading,  the / constants  are  identically  zero. 

It  is  noted  that  these  dispersion  curves  without  the  fluid  load  and  with  the  double  fluid  load 
have  both  been  previously  derived.  The  plate  dispersion  curve  without  fluid  loading  is  known  as 
the  Rayleigh-Lamb  frequency  equation  for  the  propagation  of  waves  in  a  plate,  which  is  given 


tan(/?/?/2)  I"  4  afik]  1 
tan(a/z/2)  +{(k2x  -P2f\ 


(86) 


Equation  (79)  (or  equation  (84))  with  fn  =  0  and  equation  (86),  although  not  identical,  have  the 
same  zeros  that  correspond  to  the  branches  of  the  dispersion  curves  for  the  plate  without  fluid 
loading.  The  plate  dispersion  curve  with  the  double  fluid  load  for  the  case  of  symmetrical  wave 
response6  is 


and  for  the  case  of  antisymmetrical  response6  is 
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Equation  (84)  and  equations  (87)  and  (88),  although  not  identical,  have  the  same  zeros  that 
correspond  to  the  branches  of  the  dispersion  curves  for  the  double  fluid-loaded  plate.  Figures  of 
calculated  dispersion  curves  are  shown  in  the  numerical  example  section. 
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U xs  (kx  ,z,co)  -  p3  {cos(az)  -  cos (fih)  cos[a(z  +  /?)]} 

+  p4  {cos  (/?z)  -  cos(or/z)  cos[/?(z  +  A)]} 


+ p5  sin(<2/2)sin[/?(z  +  h )]  +  p6  sin(/?/z)sin[a(z  +  h)] , 

with 

(92) 

Pi (kx , (o)  =  %ia/3kl (fi2  -k2)  , 

(93) 

Pa (kx , o>)  =  -4i ccJ3kx (fi2  -k2)2  , 

(94) 

P5(kx,6))  =  -\6ia2fi2k3x  , 

(95) 

and 

p6(kx,co)  =  2ikx(jB2  -k2f  . 

(96) 

For  the  single  fluid-loaded  plate,  the  transfer  function  of  the  normal  displacement  at  location  z 

divided  by  the  excitation  level  is  equal  to 

umm  uT;(kx,z,oj) 
PM  pA  s(kx,d)) 


(97) 


where 

UTZS  (kx  ,z,co)  =  p7  (sin(arz)  -  cos (/%)  sin[a(z  +  A)]} 

+  pg  |sin(/?z)  -  cos(ah)  sin[/?  (z  +  /?)]} 

+  pg  sin(a/2)  cos[/?(z  +  h )]  +  p10  sin (J3h)  cos [a(z  +  h)] ,  (98) 

with 

P7(kM  =  Sa2pk2(p2 - k l)  ,  (99) 

Ps(kx,co)  =  -4cck2x(/J2-k2)2  ,  (100) 

p9 (kx ,a>)  =  l6a 2  fik*  ,  (101) 

and 

pl0(kx,co)  =  2a(/32-k2xy  .  (102) 
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The  transfer  function  for  the  case  of  no  fluid  loading  is  given  in  equations  (91)  and  (97),  with  the 
/terms  set  equal  to  zero  and  each  equation  divided  by  two. 


For  the  double  fluid-loaded  plate,  the  transfer  function  of  the  tangential  displacement  at 
location  z  divided  by  the  excitation  level  is  equal  to 


Ux(kx,z,co)  UTxd  {kx,z,  a) 

Pe(o))  pAd(kx,(o) 


(103) 


where 


U™  ( kx,z,o))  =  p3  {cos(az)  -  cos(flh)  cos  [a(z  +  /*)]} 

+  p4  {cos(/?z)  -  cos  (ah)  cos  [J3(z  +  A)]} 

+  p5  sin(or/?)  sin[/?(z  +  h)]  +  p6  sin(/%)sin[a(z  +  h)} 

+  /4  sin(a&)cos[/?(z  +  A)]  +  fs  sin(/3h)cos[a(z  +  h )] , 


(104) 


f4(kx,a>)  =  4pf{yp)~' a2  fikx(p2  +  k2)2 


(105) 


f5(kx,co)  =  -2pf(ypY'akx{p2  -kx)(/32  +  k2)2 


(106) 


For  the  double  fluid-loaded  plate,  the  transfer  function  of  the  normal  displacement  at  location  z 
divided  by  the  excitation  level  is  equal  to 


Uj{kxtz,G>)  _  U?(kx,z,a>) 

PM)  /A/ (*,>«) 


(107) 


where 


=  />7{sin(az)-cos(//j)sin[a(z  +  /?)]} 

+  pt  {sin  (pz)  -  cos(ctfz)  sin[/?(z  +  A)]} 

+  p9  sin(ah)cos[j3(z  +  h)]  +  pxo  sin(//?)cos[a(z  +  h)] 
+fb  sin(a/?)sin[/?(z +  //)]  +  /7  sin(/?/z)sin[«(z  +  /7)] , 


(108) 


feiK’®)  =  -4'Pf(rpY'a2k2x(/32  +k2)2 


(109) 
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and 


f1(kx,<»)  =  -2ipf(ypr1a2(/32 -k2)(p2  +  k2)2  .  (110) 

Several  specific  transfer  functions  of  this  system  with  a  further  reduced  form  are  listed  in  the 
appendix.  These  functions  correspond  to  displacement  on  one  side  of  the  plate  divided  by  the 
displacement  on  the  other  side.  Figures  that  compare  the  transfer  functions  are  shown  in  the  next 
section  —  A  Numerical  Example. 

Once  the  closed-form  transfer  functions  are  known,  closed-form  expressions  for  the  insertion 
loss  and  echo  reduction  can  be  determined  by  inserting  equation  (108)  into  equation  (77)  for 
insertion  loss  and  into  equation  (78)  for  echo  reduction.  For  the  insertion  loss  term,  this  process 
results  in 

IL(kx,a>)  =  201og10 
where 

<t>,  ( kx  ,(o)  =  pn  (kx ,  a)  sin {ah)  +  pn  (kx ,  a)  sin(/%)  , 
with 

Pu{kx,G>)  =  *Pf<X2Pk2x(p2+k2x)2 
and 

pn(kx,co)  =  2pfa{p2  -k2x)2{p2  +k2x)2  . 

For  the  echo  reduction  term,  this  process  yields  the  expression 

CT(t„m)  =  201ogi0^(f-’^  ,  (115) 

where 

<f)E  (kx ,  co)  =  />,  (kx ,  aj)  cos  {ah)  cos  (ph) 

+  [p2  (kx  ,G))-f3  (kx ,  (o)\  sin(«/i)  sin  {ph)  -  p](kx,d)) .  (116) 
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(112) 

(113) 

(114) 


•i  rp^AKi^l 
MKM 


(ill) 


A  NUMERICAL  EXAMPLE 


A  numerical  example  is  discussed  to  illustrate  the  effects  of  fluid  loading  on  a  plate.  A 
baseline  problem  is  defined  that  corresponds  to  a  mildly  stiff  elastomeric  solid  in  contact 
with  sea  water  on  one  or  two  sides.  The  plate  material  properties  are  as  follows:  Young’s 
modulus  is  2s=108N/m2,  density  is  p  =  1200 kg/m3,  Poisson’s  ratio  is  v=  0.4  (dimension¬ 
less),  and  thickness  is  h  —  0.1  m.  The  sea  water  has  a  compressional  wave  speed  of 
cf  ~  1500  m/s  and  a  density  of  pf  =  1025  kg/m3.  The  calculated  Lame  constants  are 
^  ~  1-43  x  10  N/m2  and  p  =  3.57  x  107  N/m2 .  The  calculated  dilatational  wave  speed  is 
ca  ~  423  m/s,  and  the  calculated  shear  wave  speed  is  cs  =  173  m/s. 

Figure  4  compares  the  dispersion  curve  (equation  (79))  of  the  plate  with  a  single  fluid  load 
(solid  line)  to  the  dispersion  curve  of  a  plate  without  fluid  loading  (dashed  line).  In  figure  4 
and  the  ensuing  figures,  the  abscissa  and  ordinate  have  been  normalized  such  that 


Q  = 


ha> 


7ZC' 


and 


k  = 


(117) 


(118) 


where  Q  is  nondimensional  frequency  and  a:  is  nondimensional  wavenumber. 

Figure  5  compares  the  dispersion  curve  (equation  (84))  of  the  plate  with  a  double  fluid  load 
(solid  line)  to  the  dispersion  curve  of  a  plate  without  fluid  loading  (dashed  line).  In  general, 
most  of  the  waves  tend  to  be  slowed  down  by  the  addition  of  either  a  single  or  double  fluid  load. 
In  this  figure,  the  solid  line  originating  at  the  origin  and  terminating  at  approximately  Q  =  10  and 
at  =  1.15  corresponds  to  a  compressional  wave  in  the  fluid,*  which  corresponds  to  transfer 
function  nulls  rather  than  maximum  response. 

The  dispersion  curves  shown  in  figures  4  and  5  are  useful  for  a  number  of  reasons. 
Specifically,  they  relate  frequency  to  the  wavenumber  of  propagating  waves  in  an  unbounded 

♦In  figure  4,  this  compressional  wave  is  represented  by  the  nearly  vertical  broken  line  in  the  first  and  second 
grids  on  the  left-hand  side  of  the  illustration. 
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Normalized  Frequency  (£2  =  hm/nc 


Figure  4.  Dispersion  Curve  for  a  Single  Fluid-Loaded  Plate 
Compared  to  a  Plate  Without  Fluid  Loading 
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Normalized  Frequency  (n  =  hco/ncs ) 
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Figure  5.  Dispersion  Curve  for  a  Double  Fluid-Loaded  Plate 
Compared  to  a  Plate  Without  Fluid  Loading 


medium.  That  is,  at  any  given  frequency,  the  maximum  wavenumber  of  the  energy  traveling  in 
the  structure  can  be  determined,  which  is  advantageous  for  the  analysis  of  turbulence-loaded 
structures  where  nonacoustic  wave  energy  is  present.  Additionally,  these  curves  can  help  in 
finite  element  models  where  proper  mesh  size  is  related  to  the  wavenumber  of  energy 
propagation.  Furthermore,  these  plots  show  the  effects  of  fluid  loading  for  each  mode  of 
propagation  in  the  structure. 

Figure  6  plots  the  transfer  function  of  the  tangential  displacement  divided  by  pressure  versus 
normalized  wavenumber  with  a  normalized  frequency  of  Q  =  5  at  location  z  =  -0.04  m.* 

Figure  7  shows  the  transfer  function  of  the  normal  displacement  divided  by  pressure  versus 
normalized  wavenumber  with  a  normalized  frequency  of  Q  =  5  at  location  z  =  -0.04  m.  In  both 
figures,  the  solid  line  represents  the  single  fluid-loaded  plate  and  the  dashed  line  shows  the  plate 
without  fluid  loading. 

Figure  8  plots  the  transfer  function  of  the  tangential  displacement  divided  by  pressure  versus 
normalized  wavenumber  with  a  normalized  frequency  of  Q  =  5  at  location  z  =  -0.04  m.  Figure  9 
shows  the  transfer  function  of  the  normal  displacement  divided  by  pressure  versus  normalized 
wavenumber  with  a  normalized  frequency  of  Q  =  5  at  location  z  =  -0.04  m.  In  these  two  figures, 
the  solid  line  is  the  double  fluid-loaded  plate  and  the  dashed  line  is  the  plate  without  fluid  loading, 

Figure  10  shows  a  displacement  shape  for  the  n  =  1  antisymmetric  mode  of  the  double  fluid- 
loaded  plate,  the  single  fluid-loaded  plate,  and  the  plate  without  fluid  loading.  The  figure  on  the 
left  illustrates  plate  thickness  versus  tangential  displacement,  and  the  figure  on  the  right  shows 
plate  thickness  versus  normal  displacement.  The  solid  line  is  the  double  fluid-loaded  plate,  and 
the  dashed  line  represents  the  single  fluid-loaded  plate  and  the  plate  without  fluid  loading.  The 
displacement  shapes  were  determined  by  taking  a  point  on  the  n  =  1  branch  of  each  dispersion 
curve  for  the  separate  cases  and  then  using  these  values  to  compute  the  displacements.  For  the 
double  fluid-loaded  plate,  the  values  of  this  point  were  Q  =  2.94  and  k=  2.50;  for  the  single 
fluid-loaded  plate,  these  values  were  Q  =  3.05  and  k=  2.50;  and  for  the  plate  without  fluid 
loading,  these  values  were  Q  =  3.17  and  k=  2.50.  Note  from  equations  (91),  (97),  (103),  and 
(107)  that  the  displacement  shape  is  contained  entirely  in  the  numerator  and  that  the  location  of 


*In  figures  6  through  9,  the  top  plot  is  the  magnitude  and  the  bottom  plot  is  the  phase  angle. 
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Phase  Angle  (o)  Magnitude  (dB  ref  m/Pa) 


Normalized  Wavenumber  (xr  =  hk/n) 


Figure  7.  Transfer  Function  of  Normal  Displacement  Divided  by  Pressure 
Excitation  Versus  Normalized  Wavenumber  for  a  Single  Fluid-Loaded  Plate 
Compared  to  a  Plate  Without  Fluid  Loading 


Phase  Angle  (°)  Magnitude  (dB  ref  m/Pa) 
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Figure  8.  Transfer  Function  of  Tangential  Displacement  Divided  by  Pressure 
Excitation  Versus  Normalized  Wavenumber  for  a  Double  Fluid-Loaded  Plate 
Compared  to  a  Plate  Without  Fluid  Loading 


Phase  Angle  (o)  Magnitude  (dB  ref  m/Pa) 


Figure  9.  Transfer  Function  of  Normal  Displacement  Divided  by  Pressure 
Excitation  Versus  Normalized  Wavenumber  for  a  Double  Fluid-Loaded  Plate 
Compared  to  a  Plate  Without  Fluid  Loading 


Plate  Thickness  (m) 


Figure  10.  Displacement  Shape  of  then  =  l  Mode  for  a  Double  Fluid-Loaded  Plate 
Compared  to  a  Single  Fluid-Loaded  Plate  and  a  Plate  Without  Fluid  Loading 
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the  mode  in  the  wavenumber-frequency  plane  is  contained  entirely  in  the  denominator. 
Additionally,  because  the  single  fluid-loaded  plate  contains  no  fluid  terms  in  the  numerator,  it 
has  a  displacement  shape  identical  to  the  plate  without  the  fluid  load.  However,  the  double  fluid- 
loaded  plate  does  contain  fluid  terms  in  the  numerator,  and  thus  its  displacement  shape  is 
different  from  both  the  single  fluid-loaded  plate  and  the  plate  without  fluid  loading.  Comparison 
of  displacement  shapes  at  other  modes  yields  similar  results  to  those  of  the  unloaded,  single- 
loaded,  and  double-loaded  fluid  plate  displacement  shapes  shown  in  figure  10. 

Figure  1 1  plots  the  magnitude  of  the  insertion  loss  versus  normalized  frequency  for  the 
double  fluid-loaded  plate.  The  solid  line  corresponds  to  an  acoustic  wave  at  45°  incident  to  the 
plate,  the  dashed  line  corresponds  to  an  acoustic  wave  at  0°  incident  to  the  plate  (broadside),  and 
the  x’s  correspond  to  a  closed-form  solution  given  by20 


71(0,  <y)  =  101og10 


I 

CN 

/ - \ 

V 
1  ' 

i... 

sxa.2{kdh)  +  \ 

L  4m  J 

where 


m  = 


Pcd 


PfCf 


(119) 


(120) 


Figure  12  shows  the  magnitude  of  the  echo  reduction  versus  normalized  frequency  for  the 
double  fluid-loaded  plate.  The  solid  line  corresponds  to  an  acoustic  wave  at  45°  incident  to  the 
plate,  the  dashed  line  corresponds  to  an  acoustic  wave  at  0°  incident  to  the  plate  (broadside),  and 
the  x’s  correspond  to  a  closed-form  solution  given  by20 


££(0,<y)  =  10  log 


10 


Am2 


[(1  —  w  )  sin  (kdh)] 


•  +  1 


(121) 


It  is  noted  that  equations  (119)  and  (121)  are  valid  solely  for  excitation  at  zero  wavenumber, 
where  only  dilatational  wave  motion  is  present.  They  do  not  account  for  nonzero  wavenumbers, 
where  the  effects  of  shear  wave  motion  is  present.  Also  note  in  figures  1 1  and  12  that  the  shear 
wave  energy  originating  at  the  integer  normalized  frequencies  is  present  and  creates  a  discontin¬ 
uous  function  at  the  nonzero  wavenumber,  which  is  slightly  (positively)  shifted  with  respect  to 
the  frequency. 
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Figure  11.  Insertion  Loss  Versus  Normalized  Frequency 
for  a  Double  Fluid-Loaded  Plate 
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Figure  12.  Echo  Reduction  Versus  Normalized  Frequency 
for  a  Double  Fluid-Loaded  Plate 


CONCLUSIONS  AND  RECOMMENDATIONS 


This  report  has  derived  the  closed-form  solution  to  the  expressions  for  insertion  loss  and 
echo  reduction  using  a  sample  insonified  at  nonzero  wavenumbers.  Additionally,  the  closed- 
form  transfer  functions  were  derived  for  plate  motion  divided  by  excitation  with  fluid  loading 
on  one  and  both  sides.  The  dispersion  equation  was  also  formulated  based  on  the  matrix 
equations  and  was  then  compared  to  previously  available  dispersion  equation  forms.  It  was 
shown  that  the  zeros  of  both  dispersion  equations  were  in  the  same  location.  Furthermore,  the 
displacement  shapes  of  the  system  modes  were  determined,  and  it  was  found  that  the 
displacement  shapes  of  the  plate  modes  were  identical  for  a  plate  with  no  fluid  loading  and  for 
one  with  fluid  loading  on  a  single  side.  However,  fluid  loading  on  both  sides  produced  a 
different  displacement  shape.  Finally,  closed-form  transfer  functions  were  derived  that 
compared  displacement  across  the  plate  rather  than  referenced  to  excitation  levels.  It  was 
shown  that  plate  dispersion  curves  were  eliminated  from  the  expressions  that  resulted  when 
transfer  functions  were  derived  across  the  plate. 

Future  work  in  this  area  should  investigate  finite-sized  plates,  as  well  as  anisotropic 
response. 
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APPENDIX 


REDUCED-FORM  TRANSFER  FUNCTIONS  OF  DISPLACEMENT 

Several  reduced-form  transfer  functions  of  displacement  at  one  side  of  the  plate  divided  by 
displacement  at  the  other  side  of  die  plate  are  presented  in  this  appendix.  These  functions  are 
especially  useful  when  the  data  obtained  during  a  test  are  accelerometer  or  laser  velicometer  data 
rather  than  fluid  pressure  data.  Note  that  when  the  transfer  function  is  formulated  in  this  manner, 
the  determinant  is  divided  out  of  the  equations.  Moreover,  the  response  of  the  plate  without  fluid 
loading  is  identical  to  the  response  of  the  single  fluid-loaded  plate. 

The  transfer  function  for  the  single  fluid-loaded  plate  and  the  plate  without  fluid  loading 
corresponding  to  the  tangential  displacement  at  z  =  a  divided  by  the  tangential  displacement  at 
z  =  b  is 


Ux(kx,a,a>)  _  c, [cos(ah) - cos(fth)] 

Ux  (kx ,  b,  oj)  cos  (J3h)  cos  (ah)  +  c2  sin  (fih)  sin(a/j)  - 1 


where 


cl{kx,co) 


F+kl 

/?2-3  k2x 


(A-2) 


and 


c2(kx,o)) 


(jr-klf-ZcSpX 

2ccfi(p2  -  kx  )(/?2  -  3k x )  ‘ 


(A-3) 


The  transfer  function  for  the  single  fluid-loaded  plate  and  the  plate  without  fluid  loading 
corresponding  to  the  normal  displacement  at  z  =  a  divided  by  the  normal  displacement  at  z  =  b  is 


U z  (kx  ,a,(o)  _  sin(a/z)  +  c3  (kx ,  co)  sin (J3h) 

Uz(kx,b,co)  sm{ah)  cos(/%)  +  c3  {kx ,  o))  cos(ah)  sin(j3h) 


where 


A-l 


(A-5) 


c3(kx,co)  = 


(/?2-*2)2 

4a/3k2 


The  transfer  function  for  the  double  fluid-loaded  plate  corresponding  to  the  tangential 
displacement  at  z  =  a  divided  by  the  tangential  displacement  at  z  =  b  is 


Ux(kx,a,co)  _  Uax{kx,co) 

Ux{kx,b,G>)  Ux(kx,co)  ’ 

where 

Ux  ( kx ,  co)  =  c,  (kx ,  co)  cos  (ah)  +  dx  (kx ,  co)  sin(ah ) 

-  c,  (kx ,  co)  cos(j3h)  +  d2  (kx ,  co)  sin {ph) , 

Ux(kx,eo)  =  cos(ah)  cos  (fih)  +  d{  (kx ,  co)  sin  {ah)  cos  (J3h) 

+  d2(kx,eo)  cos  (ah)  sin  (fih)  +  c2  (kx ,  co)  sin(a/>)  sin  {fih)  - 1 , 


d](kx,a) 


~i/>/«(^2+^2)2 

rP(P2  ~k2x){p2  -3k2)  ’ 


and 


d2{kx,co) 


2ypp(p2  -3k2x)  ‘ 


(A-6) 

(A-7) 

(A-8) 

(A-9) 

(A- 10) 


The  transfer  function  for  the  double  fluid-loaded  plate  corresponding  to  the  normal 
displacement  at  z  =  a  divided  by  the  normal  displacement  at  z  =  b  is 


U:(kx,b,co)  Ubz(kx,co)  ’ 
where 

(kx  ,o))  =  c2  (kx ,  co)  sin  (fih)  +  sin  (ah)  , 


(A-l  1) 


(A- 12) 


A-2 


Ubz  (kx  ,eo)  =  c3  (kx ,  co)  cos  (ah)  sin  (J3h)  -  d3  (kx ,  co)  sin  {ah)  sin(/%) 

+  sin(a/z)  cos (J3h) ,  (A- 13) 


and 


d3(kx,aj) 


i  p,(fi2+k2x)2 
4yppk2x 


(A- 14) 


A-3  (A-4  blank) 


INITIAL  DISTRIBUTION  LIST 


Addressee  No.  of  Copies 

Office  of  Naval  Intelligence  (ONI  241  -  J.  Zilius,  T.  Morgan; 

ONI  263  -  S.  Brown  (2))  4 

Office  of  Naval  Research  (ONR  32 1  -  R.  Elswick;  ONR  333  -  K.  Ng)  2 

Defense  Technical  Information  Center  2 

Center  for  Naval  Analyses  1 


